From a nonlinear string to a weakly interacting Bose gas 
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We investigate a real scalar field whose dynamics is governed by a nonlinear wave equation. We 
show that classical description can be applied to a quantum system of many interacting bosons 
provided that some quantum ingredients are included. An universal action has to be introduced in 
order to define particle number. The value of this action should be equal to the Planck constant. 
This constrain can be imposed by removing high frequency modes from the dynamics by introducing 
a cut-off. We show that the position of the cut-off has to be carefully adjusted. Finally, we show 
the proper choice of the cut-off ensures that all low frequency eigenenmodes which are taken into 
account are macroscopically occupied. 



I. INTRODUCTION 

Description of interacting many body quantum system 
it is a very difficult task. Except of a few rather academic 
problems exact solutions are not accessible and some ap- 
proximated methods are necessary. Recent attempts of 
description of a Bose-Einstein condensate at finite tem- 
perature showed that it is possible to significantly 
simplify theoretical methods of solving quantum dynam- 
ics. This is due to the assumption of macroscopic occu- 
pation of single particle modes. It is worth noticing that 
Bose-Einstein condensate is not the only system with 
macroscopic occupation of quantum states. 

Historically the first and the best known system with 
such properties is electromagnetic field. As long as in- 
tensity is large classical approach based on Maxwell's 
equations is valid. The classical point of view was the 
only one used until the beginning of the twentieth cen- 
tury. However this approach turned out to be inadequate 
to describe experiments with small intensities and there- 
fore the concept of photon had to be introduced. On the 
other hand, physical phenomena involving macroscopi- 
cally large field amplitudes (or equivalently large num- 
ber of photons) are successfully described by the classical 
electric and magnetic fields [1| . Although quantization of 
electromagnetic field is a well established procedure, the 
inverse procedure, i.e. substitution of a quantum field by 
a classical one is often heuristic and based on physical 
intuition rather than formal arguments. At this point we 
should mention the Glauber theory of coherent states |j| 
of electromagnetic field. If the field is in a coherent state 
then electric and magnetic filed operators can be substi- 
tuted by their non-vanishing mean values in this state. 
These mean values are interpreted as classical fields. 

Although coherent states provide a link between classi- 
cal and quantum theories the situation is not that simple 
in case of particles with non-zero mass. Superselection 
rules do not allow for superposition of states with differ- 
ent number of particles. Therefore, one cannot introduce 
coherent states for such fields. Such a situation takes 
place in the case of atomic Bose condensates, where the 
number of particles is fixed. 

In Bose condensates at low temperatures only few 



lowest energy levels are macroscopically occupied. One 
might expect that discrete structure of the matter field is 
not essential and description of the system by a 'classical 
wave' should be valid. This expectation has no rigorous 
justification. However an ingenious idea of Bogoliubov 
consisting in substitution of the annihilation operator of 
a particle in the condensate mode by a c-number ampli- 
tude Q is extremely successful and widely used. This ap- 
proach leads to a mean field description of the system in 
terms of a classical fields satisfying the Gross-Pitaevskii 
equation. 

Thus the standard theory of Bose condensate at zero 
temperature is based on the Bogoliubov method. Re- 
cently this idea was extended to finite temperatures and 
is called the classical fields method. It is successful in 
describing equilibrium properties of Bose condensate at 
finite temperatures, excitations spectrum, dissipative dy- 
namics of vortices and many other finite temperature 
phenomena Q. 

The classical fields method is based on heuristic substi- 
tution of Bose operators by c-number amplitudes. This 
substitution can be easily justified at very low tempera- 
tures since practically all particles occupy the condensate 
mode. On the other hand it seems questionable at tem- 
peratures close to the critical temperature. Nevertheless 
the classical fields method works quite well up to the 
condensation point. It is known that the Bose-Einstein 
condensation is a quantum phenomenon and the critical 
temperature depends on the Planck constant. 

The Planck constant is usually introduced into the- 
ory through commutation relations for the field opera- 
tors. Those however are violated in the classical fields 
approach. Nevertheless this approach preserves some 
bosonic features of the system. Therefore it is justified 
to ask which features of the quantum system are taken 
into account in the classical field methods and which as- 
sumptions leads to correct prediction of the condensation 
temperature. In particular one may ask how does the 
Planck constant enter into the classical field method. 

In order to answer these questions, at least partially, we 
choose the following approach. We start with a classical 
system of interacting harmonic oscillators and study its 
dynamics. At a later stage we include some quantum 
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ingredients and pinpoint the moment where the Planck 
constant appears. This approach justifies the classical 
fields method and shows its limitations. 

The paper is organized as follows. In Sec. II we in- 
troduce the model. We discuss the numerical techniques 
that we use and main features of numerical solutions. A 
particular form of interactions between oscillators leads 
to nonlinear equations similar to those studied by Fermi, 
Pasta and Ulam 0|. Unlike the Fermi-Pasta-Ulam results 
the dynamics given by our model leads to thermalization 
of the system. In Sec. Ill we analyze the state of ther- 
mal equilibrium reached by the system, we study energy 
equipartition, and define the temperature of the system. 
In Sec. IV we analyze the system in terms of quasiparti- 
cles and occupation of single particle state. We show that 
elementary excitations (phonons) are distributed accord- 
ing to the low frequency part of Bose statistics. High 
frequency part is introduced by hand with the help of 
properly chosen cut-off. In this way we mimic quantum 
statistics in the whole range of frequencies. In sec. V we 
conclude by summarizing our results and showing their 
implications for the classical fields method. 

II. DYNAMICS OF A NONLINEAR STRING 

In this section we are going to introduce the model 
and its basic equations. We consider a one dimensional 
elastic string of length L and linear mass density p. In or- 
der to find dynamical equations of motion we divide the 
string into N — 1 elements of length Iq = L/N and mass 
to = lop. Each element is replaced by a point-like par- 
ticle interacting with the nearest neighbor via harmonic 
forces. The restoring force F acting on each particle is 
proportional to displacement A?o from its equilibrium po- 
sition, F = —Y(AIq/Iq) where Y is the Young modulus. 
Thus the string can be viewed as N particles moving on a 
line, each of them connected to two neighbors by a spring 
with equilibrium length lo and elastic constant K = YJIq. 
We denote equilibrium positions of each particle (oscilla- 
tor) by Xj = jlo (j — 1, ■ • • , N) and their displacements 
from equilibrium (along the axis of the string) by <f>j . The 
Newton equations of motion for the displacements are: 

mij)j = -K(2(j>j - <j> j+ i - <pj-i). (1) 

For the future convenience we assume periodic boundary 
conditions, i.e. = <j)j+N- Eq. (QJ is used in different 
areas of physics, e.g. in description of vibration of one- 
dimensional crystal lattice. Analytic solutions of Q are 
available in terms of plane waves, @. 

Let us remark, that in the limit of continuous medium, 
N — ► oo (i.e. Iq — » 0) Eq.JQ takes the form of wave 
equation: 

where c = \jYj p is the velocity of sound. In the language 
of classical field theories equation (1) describes free scalar 



field of zero charge. In the following we will use the 
discrete version of the model. 

We will now take into account a nonlinearity. For 
simplicity we assume that the nonlinear interaction is of 
short range (local) and the dynamical equation is: 

vn^i = -K(2y j - <f> j+1 - 0j-_ x ) - A$, (3) 

where A is a real parameter. This form of interaction is 
widely used in various areas of physics, in particular in 
the so called <j) 4 field theory, 0|. 

Equation 101 is very similar to the one which appears in 
the famous Fermi-Pasta-Ulam (FPU) problem Q. There 
are, however, two differences. First, in the FPU case 
displacements of the first and the last oscillators are set 
to zero (as opposed to periodic boundary conditions as- 
sumed here). Secondly, the nonlinear term in the FPU 
equation is of a different form. The authors consid- 
ered non-local nonlinear forces, for example of the form 
{4> j — <j>j-i) r , where r = 2 or r = 3. The results of FPU 
calculations show that the system shows 'very little, if 
any, tendency toward equipartition of energy among de- 
grees of freedom', Q. On the contrary, as we are going 
to show in the following, the system described by Eq.(j2ji 
reaches a state of thermal equilibrium characterized by 
equipartition of energy. 

Let us introduce natural units: 1) unit of length L, 2) 
unit of time to — Iq/c and 3) unit of energy e = KL 2 . 
The set of coupled nonlinear equations takes the form: 

4>j = -(2<£j - <p j+1 - - \$, (4) 

where 

A = ±1? (5) 

is the nonlinear coupling constant. We assume that A, 
does not depend on Iq. Therefore coefficients in Eq.(@J 
do not depend on the number of oscillators N. 
Eq. (2J leads to energy conservation: 

E = \ Y$ + - + = const - ^ 
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We solve the set of equations J3J numerically for vari- 
ous initial energies and various number of oscillators N . 
All initial displacements and velocities are gener- 
ated from uniform probability distribution on the inter- 
val [—ip,ip], where tp is a parameter. Its value has to be 
adjusted according to the value of initial energy. In gen- 
eral, the larger initial energy the larger the value of tp. 
In our numerical simulations we adjusted the time step, 
and therefore accuracy of numerical procedure, to ensure 
energy conservation. 

Due to assumed boundary conditions it is convenient 
to analyze the results of our simulations in the basis of 
plane waves 

^■(i) = -L &fc(*)e- <fci , (7) 

V k—kmin 
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where the dimensionless wave vector k takes the values 
k = nn/2 and n = —N, ...,N—1. Note that modes with 
\k\ > k max are not present in @. The cut-off of large 
|fc| is introduced in our model by discretization of space. 
We will show that the cut-off position (or equivalently the 
number of spatial grid points N) plays a very important 
role in our description. Complex amplitudes of different 
plane waves bk are subject to a constrain: 



&*(*)=*»-*(*). 



(8) 



because the displacement field <fij(t) is a real function. 
This condition is automatically satisfied in numerical im- 
plementation. 



III. THERMAL EQUILIBRIUM 

The most important observation which follows from 
numerical solutions is that the system reaches a state of 
equilibrium after some transient time. We checked this 
observation by choosing many different initial conditions 
corresponding to the same initial energy. The equilibrium 
state is characterized by randomly-looking oscillations of 
each plane wave amplitude b k {t). Its modulus |6fc(i)| 2 
oscillates in time from zero to some maximal value. The 
mode with k = has the largest amplitude. In general, 
the larger the wave vector \k\, the smaller the amplitude 
of oscillations of the corresponding mode. 

We will now analyze the stationary state. It depends 
on the energy E and the number N. Let us analyze 
the frequency spectrum of each plane wave amplitude 
b k (uj) = J dte ItA}t bk(t). It turns out that the randomly 
looking oscillations of bk(t) have a very regular spectrum. 
Typical frequency spectra of plane wave amplitudes with 
k = 0, —5, 5 are plotted in Fig. ^ 

The spectrum of each mode is composed of two peaks 
centered around uik and — uik- The peaks have finite 
width which is relatively small as compared to the cen- 
tral frequency. Therefore, in the first approximation, it 
is justified to neglect the width of peaks in the frequency 
spectrum and assume that each plane wave amplitude 
oscillates in time with two frequencies. The frequency 
ujk is defined as the weighted mean value of frequencies 
corresponding to the positive frequency peak: 



UJk 



E.>ol^HI 2 

while effective amplitude @k is: 
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OJ>0 



(9) 



(10) 



In this way we defined an unique frequency ujk and am- 
plitude pk for every mode fc. Thus the amplitude bk has 
the following time dependence: 
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FIG. 1: Spectrum of amplitudes for energy E — 0.1, A = 1, 
number of grid points N = 64 and a) k = 0, b) k = —5, c) 
k — 5. Note existence of two peaks centered at uik and — u)k- 



We checked that the frequencies ujk fulfill the following 
dispersion relation 



ujk = \/4sin - + w^, 



(12) 



where uiq is the frequency of the k = mode. Numerical 
fit shows that: 



wl=2Xa 

k=k mi 



(13) 



b k (t) = Pke 



-iiOkt 



.0* 



a iu> k t 



(11) 



and a is close to 4.5. In Fig. |2]we show the dispersion re- 
lation obtained from numerical simulations (points) and 
compare it to the analytic formula ltl2|l . Comparison 
shows remarkably good agreement. 
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FIG. 2: Dispersion relation. Points refer to numerical sim- 
ulation with E = 0.1, N = 64 and A = 1. Curve refers to 
formula il'It with a = 4.5. 
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FIG. 3: Equipartition of energy. Energy per mode, Eq. Ultifl . 
as a function of wave vector. Points refer to numerical simu- 
lation of Eq. gj with E = 0.2, N = 64 and A = 1. 



In the continuous limit, where Eq. J3J becomes a non- 
linear wave equation the 4sin 2 (fc/2) term becomes simply 
k 2 . This difference is due to discretization of space and 
the simplified form of the discrete version of the second 
spatial derivative. 

In the state of thermal equilibrium the dispalacment 
<j)j (i) can be approximated by 

1 k mal 

^(t) = -= J2 (fike~ i< - ki ~ u " ,t) +/3* fc e 4(fe -''^ fct) ). 

v N k=k mzn 

(14) 

Normal modes of the system are plane waves. Every 
mode k oscillates with two opposite frequencies LUk and 
— u)k- Amplitudes corresponding to these frequencies are 
related, \fik\ = \0-k\- Using the language of the field 
theory we can say that the positive frequency part cor- 
responds to particle-like modes while the negative fre- 
quency component corresponds to antiparticle-like exci- 
tations. 

The total energy of the interacting system can be ex- 
pressed in terms of amplitudes 0k and frequencies lou 
similarly as in the linear case: 

£ = 2 E (15) 

Factor two which appears in the above formula comes 
from the fact that each plane wave oscillates with two 
frequencies of opposite sign and each of them gives the 
same contribution to the energy. We check that disagree- 
ment between values of energy given by Eq. © and 11511 
is not larger than 5% — 10%. Therefore, the total energy 
is a sum of energies of independent modes. The hamilto- 
nian expressed in terms of plane waves is diagonal. 

One of the most important results of our numerical cal- 
culations is that total energy is evenly distributed among 
all plane wave modes. The equilibrium state reached dur- 
ing nonlinear evolution is characterized by equipartition 
of energy: 

Sk = 2uo\\(3 k \ 2 = const. (16) 



where Ek is the energy per mode. Obviously Sk = E/N. 
Fig. 10 illustrates the equipartition of energy. Let us re- 
mind that no equipartition of energy is observed in FPU 
problem. 

Having the equipartition of energy we can define a tem- 
perature of the system: 

f = s k , (17) 

where temperature is expressed in unit of e/feg (&b is the 
Boltzmann constant). 



IV. PHONONS 

In the previous section we have shown that the non- 
linear hamiltonian © can be expressed in the diagonal 
form l(T5)l for each particular state of thermal equilibrium. 
Note that it cannot be done 'globally', i.e. independently 
of the total energy of the system. This is because the 
eigenfrequencies depend on the energy. 

The diagonal hamiltonian 111 oil is a sum of energies of N 
independent harmonic oscillators of frequencies u>k • The 
interaction strength A enters this hamiltonian through 
the eigenfrequencies cuk only. 

This picture is purely classical. In what follows we will 
reformulate our results in the language of quantum many 
body theory. Our goal is to reach a better understanding 
of foundations of the classical fields method. 

The hamiltonian 11511 can be quantized with the help 
of canonical quantization procedure, i.e. by defining po- 
sition and momentum in terms of amplitudes /3k, and 
then imposing canonical commutation relation between 
them. However, bosonic commutation relations are vio- 
lated in the classical fields method - classical amplitudes 
do commute. Thus we will choose a different approach. 

First we will define dimensionless amplitudes. These 
amplitudes are the only dynamical objects appearing in 
the classical field method. In order to define them we 
have to rewrite Eci JlfiJl using quantities with their proper 
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dimensions: 



n 



n\\A k \ 



(18) 



k=k„ 



where flk = (1/ioVfc is the frequency and A k — V2L(3 k 
is a classical amplitude of each harmonic oscillator. To 
define dimensionless amplitudes we need some universal 
quantity which has dimension of an action. Such a quan- 
tity does not exist in the classical theory. Therefore we 
have to introduce this constant by hand. We will use 
symbol h for this elementary action. At this moment its 
value can be arbitrary. So far we have not mentioned 
any physical condition which could, at least in principle, 
determine its value. At a later stage this constant will 
be identified with the Planck constant. 
Let us define units of amplitudes: 



MQ = 



mQ k 



(19) 



and dimensionless amplitudes for each oscillator: 



W2/3* _ /2mtt fe 

a h — , Ti \ — \ ^ L'Pki 



A (k) V h 

and express the Hamiltonian Eci. ill8|) in terms of B k : 



(20) 



(21) 



This form of the Hamiltonian En J21ll is a familiar one. 
The energy of each plane wave mode is equal to the 
energy of elementary quantum of a given frequency Q, k 
times some positive real number: 



N k = \B k \ 



(22) 



which can be called a number of phonons. We want to 
stress once more that number of phonons defined above is 
not an integer. Therefore it can have a physical meaning 
of the number of quasiparticles only if its value is large 
as compared to one. Moreover, the value of N k depends 
still on the (arbitrary) value of the elementary action H. 
Because there is no limitation on the maximal amplitude 
of a harmonic oscillator, the value of N k can be arbitrar- 
ily large. Our classical field cannot therefore correspond 
to particles obeying fermionic statistics. However it can 
correspond to highly excited Bose field. 

Equipartition of energy Eqs. iflfil) . ifTTjl . discussed in 
the previous section, can be expressed in terms of N k : 



hfl k N k = k B T 7 



(23) 



where kg is the Boltzmann constant and fc^T = eT 
where T is known from our numerical simulations. On 
the other hands if N k corresponds to the number of 



phonons (in the 'classical limit') than in the thermal equi- 
librium it should obey Bose statistics: 



N k = 



1 



e h{Q. k -n)/k t 



1 



(24) 



where /i is a chemical potential. In the limit of low fre- 
quencies Cl k> the equilibrium occupation of the mode k 
can be approximated by: 



N k h(n k -fi) = k B T. 



(25) 



Comparison of Eq. 12511 and Eq. i'2'M) shows that occupa- 
tion of normal modes obtained in our calculations agrees 
with a low frequency limit of Bose statistics provided that 
fi = 0. 

Equipartition of energy defines the energy per mode. 
All modes with wave vectors |fc| < k max are occupied by 
phonons while other modes are empty because of the mo- 
mentum cut-off used in the implementation of the model. 
This distribution of energy is a ID analogue of the Planck 
distribution of the blackbody radiation. The Planck dis- 
tribution says that energy density grows with frequency 
up to some maximal value and then falls exponentially 
to zero. This initial growth is related to the phase space 
volume which increases with k as cx k 2 and is absent in 
the ID model studied here. The exponential decay of the 
blackbody energy density in replaced by a sharp cut-off 
at k — k max in our calculations. 

We will use the similarities described above to deter- 
mine the value of h which will allow us to find the abso- 
lute value of the number of phonons N k . Following the 
Planck idea we determine the value of h by equating the 
temperature of the system to the energy corresponding to 
the position of the maximum in the energy distribution. 
In our case this 'maximum' corresponds to the cut-off: 



= k R T. 



(26) 



Eq. l(26j) plays a crucial role in establishing a link be- 
tween classical fields approach and the quantum theory 
of Bose system. It gives the value of the Planck con- 
stant up to a numerical factor of the order of one. This 
fact establishes a limitation of accuracy of the method, 
in particular the accuracy of the number of particles, or 
equivalently the value of the temperature of the system. 
By comparing Ea J26|l and Eci. l|23ll we see that number 
of particles occupying the mode with the largest energy 
is equal to one: 



Nk 



1. 



(27) 



This means that occupation of all modes with k smaller 
than k max is macroscopic, N k > 1, which is in agreement 
with a common understanding of the classical limit of the 
quantum field. Indeed in the classical fields method the 
number of modes used in numerical implementation has 
to be carefully adjusted. As suggested in ^3] it should 
be such that the occupation of the highest mode is equal 
to 1. 
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FIG. 4: Relative number of phonons in k = mode as a 
function of temperature for different total number of quasi- 
particles: Nph = 85 (gray circles) and Nph = 175 (black 
circles). Points refer to numerical simulations of Eq.(4) with 
A = 1. Lines represent a fit with a = 4.5. 



In what follows we give an illustrative example of ap- 
plication of the defined above procedure of 'quantization 
of a scalar field' studied in this paper. We solve the dy- 
namical equations for various total energy of the system 
and various number of grid points N. In our calculations 
we used rather small values of N ranging from N = 16 up 
to N = 1024. We check if the system reached the state 
of thermal equilibrium by determining a frequency spec- 
trum ojk and amplitudes (3k at various stages of time evo- 
lution. When the stationary state is reached we control 
equipartition of energy and determine the temperture. 
Finally, occupation of different eigenmodes is obtained 
from the following relation: 



\0k\ 



(28) 



V. CONCLUDING REMARKS 

In this paper we started from the description of 
a purely classical system. Our numerical simulations 
proved that nonlinear interactions drive the system to- 
wards an equilibrium. Energy in the equilibrium is evenly 
distributed among all modes of \k\ < k max . Modes cor- 
responding to |fe| > k max are not excited and do not 
contribute to the total energy. They are absent in our ap- 
proach because of the cut-off which is an essential ingre- 
dient of the method. In order to 'quantize' the dynamics 
we introduced the Planck constant. Then we could de- 
fine dimensionless amplitudes of each oscillator and the 
number of excitation quanta of each mode - phonons. 
However, the number of phonons depends on the value 
of h which, at this stage of the approach, can be arbi- 
trary. In order to assign a value to h we followed the 
Planck's idea. We used the fact that Planck's distribu- 
tion of energy leads to its equipartition for small values 
of k and exponential drop for large k. Energy distribu- 
tion reaches a maximum at frequency u>k ma x close to this 
satisfying the relation Swfc max = ksT. In our approach 
the same condition was used to determine the value of h 
and thus the number of phonons Nph . 

We believe that our studies shed more light on the 
classical fields method used for description of a weakly 
interacting Bose-Einstein condensate. The way of reason- 
ing goes the opposite way in 'derivation' of this method. 
One starts with quantum many body theory of a Bose 
gas with short range interactions. Two particle interac- 
tion energy is g/V where g is proportional to the s-wave 
scattering length and V is the volume. The number of 
macroscopically occupied modes is being assumed a pri- 
ori by choosing a value of the cut-off momentum k max . 
Annihilation operators of these selected modes are then 
substituted by classical amplitudes ak and Heisenberg 
operator equations transform into equations for classical 
amplitudes ak'. 



Evidently Nk ma:c — 1, what is equivalent to the condi- 
tion (l26l) . If total number of phonons Nph = Nk is 
different than assumed we have to repeat the calcula- 
tions changing the number of grid points. In general the 
smaller the energy the smaller N has to be used in order 
to keep the same number of phonons. This way we are 
able to get occupations of all normal modes as functions 
of temperature for a fixed total number of phonons. 

In Fig.@| we present occupation of k = mode 
as a function of temperature for nonlinear interaction 
strength A = 1 and two different total number of 
phonons. The solid line represents a fit to the numer- 
ical results. The results show that relative occupation of 
the spatially uniform k = mode grows quite rapidly at 
temperatures close to zero. There is no phase transition 
in the studied system because, first of all, our system 
is one-dimensional, and also elementary excitations are 
massless. The main goal of these calculations was to il- 
lustrate how the classical field method works in practice. 



in—ctk 
dt 



2 1.2 



2m 



gNp 
V 



22 ot*ki a ^ a ki-k 2 +k, (29) 
fci ,fe 2 



where N p is the total number of particles. This quan- 
tity corresponds to the number Nph of phonons i the 
discussed classical model. Note that N p does not enter 
dynamical equations alone - it is multiplied by the inter- 
action strength g. 

Relative occupation of different k modes is closely re- 
lated to the classical amplitudes: 



Nk 



(30) 



Let us observe that the Planck constant appears in the 
equations right from the beginning. But its value can 
be arbitrary - the classical amplitudes ak are a direct 
analogue of amplitudes Bk of our model. 

Dynamics of classical fields leads to the equipartition 
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of energy: 

huJknk — ^ — const. (31) 
l\p 

Only the value of hw^n^ is known from numerical cal- 
culations, therefore ll-'illi allows for determination of the 
ratio of ksT/N p but not the values of T, N p separately. 
Note that both N p and g are not uniquely determined be- 
cause only product gN p is an initial control parameter. 
According to the present studies this problem can be re- 
solved if the value of h is determined by the requirement 
that the energy distribution of classical fields mimics the 
quantum distribution - the position of maximum ought 
to be approximately given by 

hu kmax = k B T. (32) 

This equation together with the equipartition relation 
allows for determination of the particle number: 

N p = — !— . (33) 

In the classical fields method the above condition is jus- 
tified on a basis of heuristic arguments by saying that 
all classical modes have to be macroscopically occupied. 



The approach we used gives a new interpretation of this 
reasoning. 

Moreover, our studies show limits of the 'predictive 
power' of the classical fields method. Because the Planck 
constant is determined with limited accuracy all predic- 
tions about the total number of particles or temperature 
of the system are accurate up to a numerical factor of the 
order of one. It seems therefore that attempts of more 
accurate determination of the temperature [ll| cannot be 
free of the ambiguity related to the position of the cut- 
off. In addition the classical fields method cannot be used 
for description of very subtle effects such as a shift of the 
critical temperature of Bose-Einstein condensation due to 
interactions. For a typical system this shift is very small. 
Some efforts Q of obtaining the shift of critical tem- 
perature using classical fields method must be inevitably 
biased by the approximate character of Eo . (l32l) . 
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